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Scope 

This is the Annual Report #3 for Grant NAGS- 1605 entitled “The Propulsive Small 
Expendable Deployer System (ProSEDS)” prepared by the Smithsonian Astrophysical 
Observatory for NASA Marshall Space Flight Center. This report covers the period of 
activity from 1 August 2001 through 30 July 2002. The technical officer for this grant is 
Leslie Curtis at NASA MSFC. 
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Summary 

This Annual Report covers the following main topics: 

1. Updated Reference Mission 

The reference ProSEDS mission is evaluated for an updated launch date in the Summer 
of 2002 and for the new 80-s current operating cycle. Simulations are run for nominal solar 
activity condition at the time of launch and for extreme conditions of dynamic forcing. 
Simulations include the dynamics of the system, the electrodynamics of the bare tether, the 
neutral atmosphere and the thermal response of the tether. 

2. Evaluation of power delivered by the tether system 

The power delivered by the tethered system during the battery charging mode is 
computed under the assumption of minimum solar activity for the new launch date. 

3. Updated Deployment Control Profiles and Simulations 

A number of new deployment profiles were derived based on the latest results of the 
deployment ground tests. The flight profile is then derived based on the friction 
characteristics obtained from the deployment tests of the F-l tether. 

4. Analysis/estimation of deployment flight data 

A process was developed to estimate the deployment trajectory of the endmass with 
respect to the Delta and the final libration amplitude from the data of the deployer turn 
counters. This software was tested successfully during the ProSEDS mission simulation at 
MSFC EDAC. 
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1. Updated Reference Mission 
1.1 Introduction 

We recall in this brief introduction a few general information about the Propulsive Small 
Expendable Deployment System (ProSEDS) in order to spare the reader from looking up 
this basic information in other documents. ProSEDS will carry out a demonstration of a 
bare electrodynamic tether for propulsion. The system will fly as a secondary payload on a 
Delta II and it will be deployed from the 2 nd stage (see Fig. 1). The electrodynamic forces 
generated by the current flowing in the conductive tether are expected to strongly increase 
the decay rate of the Delta stage. The reader should consult references 1 2 3 4 5 for a more 
detailed description of ProSEDS and the principles of operation of bare-tether anodes. 



e 


Figure 1 Schematic of ProSEDS on Delta 2 nd stage 
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Electrons are collected along the upper portion of the bare tether which is positively 
biased with respect to the surrounding plasma. They are then reemitted into the ionosphere 
at the Delta stage by an active hollow cathode device. The (electronic) current flowing along 
the tether from its top to the bottom interacts with the Earth’s magnetic field to produce an 
electrodynamic drag force which decreases the orbit of the Delta stage. 

The performance of ProSEDS will be assessed on the basis of the decay rate of the 
Delta stage which is affected (for a given tether design) by the plasma conditions at mission 
time. The launch date of ProSEDS has changed throughout its development with the latest 
launch date being in the Summer of 2002. Moreover, the modes of operation of the tether 
current (i.e., the operating cycle) have also changed, as explained later on, resulting in a 
strongly increased decay rate. 

1.2 New Mission Parameters and Current Operating Cycles 

The ProSEDS mission had some notable changes in terms of launch date, starting 
altitude and current operating cycle. The launch date was moved to the Summer of 2002, 
the starting altitude was decreased from 400 km to roughly 360 km and the secondary 
operating cycle (with 80-s period) now consists of 50-s battery charging and 30-s shunt 
mode (see Fig. 3). The original primary cycle (with 60-s period as shown in Fig. 2) is 
utilized during the first 5 orbits while the secondary cycle is utilized thereafter. 

The reference mission for the scheduled launch date of July 25 lh , 2002 is analyzed in the 
following. The new launch date impacts the system performance because the current in the 
tether is a function of the plasma density which, in turn, is a function of the solar activity. 
Consequently, the electrodynamic drag force and the reentry time change because of the 
changed plasma conditions. We are presently in solar cycle 23 with the solar activity going 
down towards a minimum of activity currently estimated to happen around the years 2007- 
2008. More importantly the ProSEDS tether current operating cycle has been changed. 
The older cycle consisted of 7 orbits of primary operating cycle followed by a secondary 
cycle. The old secondary cycle consisted of 35-s open circuit mode, 5-s resistor mode, 5-s 
shunt mode and 35-s battery charging mode. The new secondary operating cycle consists 
of 30-s shunt mode and 50-s battery charging mode. The change has a positive effect (as 
shown later on) on the decay rate of the orbit because there are no periods in the secondary 
operating cycle in which the tether current is off (open circuit mode) or severely limited 
(high resistor mode). 
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Figure 2 Operating cycle #1 (60-sec cycle) - Operation on Primary Battery 
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Figure 3 New operating cycle #2 (80-sec cycle) - Operation on Secondary Battery 


Mode la is open circuit, plasma contactor OFF 

Mode lb is open circuit, plasma contactor ON 

Mode 2 is SHUNT mode 

Mode 3 is RESISTOR mode 

Mode 4 is BATTERY CHARGE mode 

Operating cycle #1 is on primary battery (first 5 orbits) 

New operating cycle #2 is on secondary battery (after the first 5 orbits) 

In the following, we show the results of the new reference mission simulation. The 
initial conditions and orbital parameters assumed for the reference simulation are as follows: 

Launch date and time: 

25 July 2002; UTC = 00: 15:00 hr:min:s (at first perigee pass) 

Orbital elements (from Boeing): 

Semimajor axis = 6736.556794 km 
Eccentricity = 0.001305 

(Initial orbit approx: 350 km X 367 km, with respect to a spherical Earth) 

Inclination = 35.363 deg 
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Right Ascension (RAAN) = 60.296 deg 
Argument of perigee = 289.996 deg 

Simulation started at first perigee pass, i.e., true anomaly = 0 
Environmental conditions: 

Nominal solar activity (50% percentile): Rz = 94 (SunSpot Number); IG = 132 
(Ionospheric Global Index); F10.7 = 141 (radio flux), Ap = 15 (geomagnetic 
index). 

System parameters: 

Satellite mass = 21.4 kg; 

Delta mass = 994 kg 

Tether mass = 10-kg wire and 1.5-kg Dyneema 
Wire external diameter = 1.2 mm 

Dyneema average cross section: ~ 1.2mm X 0.2mm 
Tether optical properties: 

Dyneema: a s = 0.1, e IR = 0.5; 

C-COR coated wire: a s = 0.9, = 0.8. 

Tether mechanical properties: EA = 15,000 N; E’A ~ 2000 Ns. 

Wire ohmic resistance = 250 ohm at 20 °C 

Electrical operating cycle: 

The primary cycle is activated at Delta Time = 12,666 s which corresponds to 4700 
s after simulation start (when ProSEDS is off the coast of Brazil). The secondary 
cycle is activated after 5 orbits, i.e., at about 32,200 s after simulation start. 

Results of the new reference mission simulation are shown in the following figures. 
Although these results have been computed for July 25, 2002, they are also valid throughout 
the Summer of 2002 because the ionospheric indexes change only a little over those months 
as explained in the next subsection. 
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ProSEDS 250 ohm@20 C, 350x367km, nom. solar, UTC 00:15:00, 25 July 2002 



Figure 4 Simulation results for Summer 2002 launch; nominal solar condition 
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ProSEDS 250 ohm@20 C, 350x367km, nom. solar, UTC = 00:15:00, 25 July 2002 
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Figure 5 Simulation results for Summer 2002 launch; nominal solar condition 
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ProSEDS 250 ohm@20 C, 350x367km, nom. solar, UTC = 00:15:00, 7/25/2002 
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ProSEDS 250 ohm@20 C, 350x367km, nom. solar, UTC 00:15:00, 25 July 2002 



Figure 7 Simulation results for Summer 2002 launch; nominal solar condition 
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ProSEDS 250 ohm @20 C, 350x367km, nom. solar, UTC 00:15:00, 25 July 2002 



Figure 8 Simulation results for Summer 2002 launch; nominal solar condition 
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The results show that the new secondary operating cycle almost doubles the decay rate 
of ProSEDS with respect to previous scenarios. The system dynamics is also more 
strongly excited with the excitation increasing steeply after one day in orbit. The excitation 
is the result of a complex skip-rope motion in which the in-plane and out-of-plane motion 
are coupled as was studied in Refs. 6 . The most noticeable indicator of this motion is the 
distance from the Delta to the end mass which is plotted in Figure 6. The fact that this 
distance does not go to zero suggests that the motion is not a planar motion but rather a 
coning motion resulting from a coupling between the in-plane and out-of plane degrees of 
freedom. The system dynamics is quite strong as it can be seen from Figure 6. The tension 
spikes in the tension plot result from recoiling after the tether has lost tension while there 
was a substantial shortening of the distance between the end mass and the Delta. The 
system, however, remains stable and it never topples over. The out-of-plane libration 
amplitude builds up to somewhat less than 30° at around 5.5 days and then decreases again; 
most likely because of an energy transfer between degrees of freedom. 

If the tether were to remain intact (neglecting the effect of atomic oxygen and possibly a 
fatal micrometeoroid impact) the Delta stage would reenter the atmosphere before 8 days. 
However, this situation is only theoretical because the atomic oxygen will chew the 
Dyneema tether before the system encounters the dense atmosphere. The area below 250 
km is were the Dyneema degradation due to atomic oxygen will become strong. The 
Dyneema tether will likely fail below that altitude, leaving the electrodynamic system is an 
uncertain state of functionality. Consequently, the rate of decay will slow down 
significantly after the tether severance with respect to the simulated decay rate that does not 
take a possible tether cut into account. 


1.3 Applicability of results to later launch dates 

We are presently in the solar cycle 23 during which the solar activity peaked in April- 
June 2000. Consequently, the solar activity and the plasma density (that is a function of the 
solar activity) will decrease over the next few years. The plasma density in the International 
Reference Ionosphere (IRI95) model is a function of two solar indexes (other than the point 
location): the Sun spot number (Rz) and the Ionospheric global index (IG). Figure 9 shows 
the prediction of the Sun spot number for this solar cycle and the next inclusive of high 
(95% percentile) and low (5% percentile) predicted values. Figure 10 shows the only 
available prediction of the Ionospheric global index (IG). 


15 



SAO Annual Report #3, NASA Grant NAG8-1605 



Figure 9 Sunspot number predictions [NASA MSFC] 



Figure 10 Ionospheric Global Index (IG) prediction [Rutherford Appleton laboratory] 

The change of these two indexes over the 3 months of the Summer 2002 is fairly 
moderate and consequently the result of the July 25 th simulation are valid in good 
approximation throughout the Summer of 2002. 
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1.4 Extreme electrodynamic forcing 

For the purpose of evaluating the system behavior and stability under maximum forcing 
of the electrodynamic drag a case was run for a plasma density that is (artificially) twice the 
plasma density under nominal conditions. 

The double plasma density increases the decay rate to 23 km/day averaged over six days 
of secondary cycle operation. Figures 11-15 show the response of the system under these 
extreme conditions of plasma density. It is notable that the decay rate does not double with 
respect to the baseline case thanks to the ability of the bare tether to adjust in part to 
changing plasma conditions. 

The dynamic response of the system is stable. The growth of the in-plane libration is 
actually less than for nominal solar conditions. The higher growth of the in-plane libration 
in the nominal case was likely due to a transfer of energy from the out-of-plane to the in- 
plane degrees of freedom as it can be inferred from Figure 6. This points again to the 
strongly non-linear characteristics of this system dynamics. The loss of tether tension most 
likely plays a key role by changing the phase relationship among degrees of freedom and 
forcing terms after the tether rebound and consequently transferring energy from some 
oscillations to others in a fairly random fashion. 

The conclusion is that the overall system stability is not impaired by a doubling of the 
plasma density. 
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Figure 12 Results for Summer 2002 launch, plasma density twice nominal 
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ProSEDS 250 ohm@20 C, 360km, double plasma, 00:015:OOT, 25 July 2002 



Figure 14 Results for Summer 2002 launch, plasma density twice nominal 
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ProSEDS 250 ohm@20 C, 360km, double plasma, 00:15:00UT, 25 July 2002 



Figure 15 Results for Summer 2002 launch, plasma density twice nominal 
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1.5 Concluding Remarks 

The changes in the operating cycles and most notably the elimination of the open-circuit 
mode during the 80-s cycles have changed the decay rate and dynamic response of 
ProSEDS. The present estimates of the orbital decay rate for a launch on July 25 th , 2002 
are 11 km/day while operating the 60-s cycle (for 5 orbits) and 19 km/day while operating 
the 80-s cycle during the remainder of the mission but before the neutral drag overpowers 
the electrodynamic drag. The decay rates estimates are valid with good approximation for 
any launch date during the Summer of 2002 because the estimated ionospheric and solar 
indexes that drive the plasma density vary only moderately over the three-month Summer 
period. 

An extreme case with a plasma density twice as high as for nominal solar conditions 
was run to check the dynamical stability of the system. The decay rate increases only by 
20% with respect to the nominal case thanks to the characteristic of the bare-wire anode to 
adjust to variations of the plasma density. The dynamics of the system remains stable. 
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2.0 EVALUATION OF POWER DELIVERED BY THE TETHER 

2.1 Introduction 

The elimination of the open-circuit mode in the 80-s operating cycle changes the amount 
of electrical energy available to recharge the battery. As a matter of fact, the change in the 
operating cycle increases the time spent in battery charging mode (i.e., mode 4) from the 
previous 35 s to 50 s per cycle. Moreover, the value of the tether resistance of the flight 
tether is lower (250 ohm at 20 °C) than the value of 265 ohm (at 20 C) utilized in the 
previous estimates of the current produced by the tether. Consequently, a new estimate of 
the tether current during battery charging was carried out as shown in the following 
subsections. 

2.2 Numerical Results 

The estimate of the tether current produced during battery charging was computed for 
minimum solar activity conditions predicted for July 2002. The minimum solar conditions, 
corresponding to a 5 % percentile probability, provide a conservative estimate of the tether 
current (and consequently power) available to recharge the secondary batteries. All other 
orbital and system parameters are the same of the reference mission described in Section 
1.0 of this report. Figure 16 shows the instantaneous tether current during battery charging 
(mode 4) and the 12-hour running average of the current. The average is computed over the 
mode-4 period of 50 s and not over the overall cycle of 80 s. 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 


Time (s) 

Fig. 16 Tether current during battery charging mode: instant value and 12-hour average 
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The average value of the current is 1.33 Amp during the first day while the peak value of 
the current is 3.7 1 Amp. The corresponding values for nominal solar conditions (the plot is 
not shown) are 1.42 Amp for the average and 3.85 Amp for the peak value during the first 
day. The average value of the current decreases by roughly 30% after 1.5 day mostly due to 
tether dynamics. 

Although the average values of the current during battery charging are similar to those 
computed for earlier launch dates, the net advantage of the new 80-s cycle is that the 
electrical energy provided by the tether (during normal battery charging operation) is higher 
because mode -4 is now 50-s long as opposed to 35 s for the older version of the cycle. 


2.3 Concluding Remarks 

The change in the 80-s operating cycle has improved the power budget for ProSEDS 
because of the longer time spent in battery-charging mode (mode 4). The estimate of the 
average current for the first day of operation in battery charging mode is 1.33 Amp for 
minimum solar conditions with the average computed over the 50-s battery charging mode. 
The estimated average value for nominal solar conditions is 1.42 Amp. The average values 
decrease by roughly 30% after 1.5 day due to system dynamics. 

The value of the average current is similar to what was computed for the previous launch 
dates. However, the total electrical energy available for battery charging (during normal 
battery operation) is better than with the older 80-s cycle because the battery charging mode 
is now 50-s long as opposed to 35 s. 
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3.0 Updated deployment control Profiles and Simulations 

3.1 Introduction 

We recall here that the ProSEDS control law consists of three distinct modes of 
operations which are activated during the deployment of the three different tether sections. 
The non-conductive 10-km-long Dyneema tether is deployed according to the SEDS-II 
feedback-feed-forward control law. During the deployment of the 4.85-km conductive wire, 
the brake is simply kept at a constant, low value in order to limit the deployment velocity. 
During the deployment of the 215-m insulated tether section, the brake is commanded to 
follow a time-based profile to slow down the deployment velocity at the end of the tether. 

As explained in details in the Annual Report #1, the control law utilizes a set of control 
parameters and a reference table that provides the feed forward information to the first 
portion of the control law. The reference table provides the control law with the nominal 
deployed length, speed (in terms of turn counts and turn count rate) and brake profiles that 
the system should follow under ideal conditions. The feedback, then, adjust the nominal fed 
forwarded brake profile based on the errors of the actual length and speed with respect to 
the nominal length and speed profiles. 

The second and third portions (for the CCOR wire and insulated section deployment) 
are open-loop control. The second portion is a tension-offset control in which the brake is 
kept constant at a low value of turns simply for increasing the tether tension and reduce the 
maximum exit velocity. The wire is coated with a fairly delicate coating that can be rubbed 
off by excessive friction. Consequently, it is not possible to utilize a feedback control that 
ramps the brake up and down. The offset value of the brake utilized during this portion is a 
constant 0.5 brake turn for the flight control law. 

In the third portion of the control law, the brake is made to follow an open loop 
rampup-constant-ramp down profile. This control law acts on the 215-m-long section of 
the tether with a quick ramp up phase. Because of the absence of a sensor that measure the 
exit velocity directly in the SEDS hardware, the velocity must be computed numerically 
from a noisy signal and then filtered to make it usable for a feedback control law. This 
process is actually used in the first portion of the control law where a delay in the 
computation of the velocity does not affect the performance of the controller. Due to the 
short reaction time available for the tether slow down, we opted for an open-loop control 
during this phase and we shaped the slow down profile in such a way that it is rather 
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tolerant of changes in the friction characteristics of the tether (see Annual Report #1 for 
more details). 

3.2 Friction model 

We rewrite in the following the frictional tension model (derived originally by J. 
Carroll, see also Ref. 7 ) of the tether and its parameters for the reader’s convenience: 

T = (r 0 + / • p ■ L 2 ■ A*,) ■ k 2 ^ e ^ e • e /l e o-0| (1) 

In equation (1), the term in round parenthesis is the frictional model of the 
tether/deployer. The first exponential function with base k is the model of the brake and the 
second exponential function with base e is the model of the tether exit guide. The model 
parameters are: 

A = 1 

rel 

A sol = annulus solidity of tether 

L = length of tether deployed 

L fin = final length of tether 

B = 2nfn (n is the number of tether turns wrapped around the brake post) 

T 0 = minimum tension 

P = linear density of tether 

I = inertia multiplier 

L = tether exit speed 

p = tether’s exit angle with respect to the local vertical 

= tether deployment null angle (orientation of the longitudinal axis of 

deployer with respect to the local vertical) 
f = friction coefficient 

n = number of brake turns 

effe = brake effectiveness coefficient 

k = base of the brake power law 

The updated values of the friction parameters obtained from the deployment tests on 
ground of the F- 1 tether are as follows: 
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Date: 8 January 2002 

From: E.C. Lorenz ini (SAO) 

Friction Parameters of the Fl flight tether 

values of the tension parameters for the 
ProSEDS tension model are as follows. These values are based 
on the deployment test (room temp.) of the Fl tether. 

Flat Dyneema (cleaned) 

T . = spectra minimum tension 

min ^ 

den = spectra linear density 

I = inertia multiplier 

f = spectra friction coeff 

E = area exponent 

ef f e 

effe = brake effectiveness Turn 
AnSol= annulus solidity 
k = base of braking power law 


= 18 mN 
= 0.2 kg/km 
= 2.0 
= 0.17 
= -0.4 
= 1.2 

= 0.2 (for L f = 10km) 
= 1.7 


Wire 

(CCOR) 


T • 

min 

= wire minimum tension 

= 75 mN 

den 

= wire linear density 

= 1.935 kg/km 

I 

= inertia multiplier 

= 1.7 

f 

= friction coeff 

II 

o 

\D 

E 

= area exponent 

= -0.65 

effe 

effe 

= brake effectiveness Turn 

= 1.2 

AnSol 

.= annulus solidity = 0.947 

(from 1.25km -> 

k 

= base of braking power law 

= 2.72 (= e) 


6.25km) 
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Insulated (Kevlar overbraided) 


T roin = insulated minimum tension 

= 

300 mN 

den = insulated linear density 

- 

3.176 kg/km 

f = friction coeff 

= 

0.22 

I = inertia multiplier 


1.7 

E = area exponent 

= 

-0.65 

effe = brake effectiveness Turn effe 


0.9 

AnSol= annulus solidity 

= 

0.947 

k = base of braking power law 

- 

2.72 (= e) 


Subsequent tests of the F-l run in hot and cold conditions validated most of the values 
of the friction parameters with the exception of the minimum tether tension which is more 
sensitive to temperature variations and even more importantly to the amount of residual 
sizing agent in the non-conductive tether (level of cleanliness). The flight tethers have all 
been cleaned in order to provide values of the minimum tension below 20 mN. However, 
the control law is tested (as shown later on) also for much higher values of the minimum 
tension. The typical range of variation of the non-conductive tether minimum tension is 
assumed conservatively to be in the range 5mN-30mN which are the limit values measured 
during the deployment tests of the many tethers developed for ProSEDS (and not just the 
F-l). 


The friction parameters of the whole tether (with the three different sections) are utilized 
to derive the reference table, that is, the reference deployment profile and a brake profile for 
the whole tether. The brake actuation is then adjusted by the feedback control law during 
the deployment of 10-km Dyneema portion while the reference brake profile is followed 
(without adjustments) during the wire and insulated portions of the tether. 

3.3 Flight reference table 

Several new reference tables were derived for the deployment control law as tether 
friction data became available during the year from the deployment tests and spooling data 
from Tether Applications, Inc. We will not describe here all the new updates but rather 
concentrate on the flight reference profile (Reference #78). 
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The process for deriving the reference profile is an optimization process that aims at 
minimizing a cost function which encapsulates the departure of the final dynamic state from 
the desired state. The desired final state at the end of deployment is for the system to be 
aligned and swing less with respect to the local vertical with a residual longitudinal velocity 
greater than 3 m/s before the beginning of the insulated portion of the wire (last 215 m of 
tether). The residual velocity is then reduced by a final activation of the brake immediately 
after the exiting of the insulated wire is sensed. The process for deriving the reference 
profile is described in more detail in the Annual Report #2. 

The flight reference table was derived based on the friction data of the F-l tether 
deployment test at room temperature and for a minimum tether tension equal to 20 mN. 
This minimum tension is slightly higher than the 18 mN measured in that deployment test 
but the 20 mN value provided a reference profile with better features (e.g., smoothness of 
the deployment velocity and progressive ramping up of the brake) than what could be 
obtained with the 18 mN value. The friction data from the room temperature test were 
preferred to the hot and cold test data because the deployer will be close to room 
temperature at the start of deployment. 

The most recent parameters adopted for deriving the flight reference profiles are as 
follows 

Orbital and ejection parameters 
Orbit: 360x360 km 
Orbital inclination: 36 deg 
Ejection velocity = 2.74 m/s 

Ejection angle = 5 deg (forward of LV with an upward deployment) 

System parameters 
Satellite mass = 21.4 kg 
Delta-II Mass = 994 kg 

Tether lengths: 10132 m Dyneema, 4865 m CCOR and 215-m insulated 

Tether length is then converted to tether turns (and velocity to turn rate) which are the 
variable utilized by the control law. Consequently, two tethers appear very similar to the 
control law if the number of turns of the three segments are about the same. The following 
files (supplied by J. Carroll, Tether Applications, Inc.) provide the spooling data for the 
F-2 and F-3 tethers, respectively 
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F-2 flight tether spooling 


Tether Applications, Inc. 

Here are the quadratic equations for the Kevlar+Dyneema , CCOR, and 
overbraided insulated wire. The values are derived from the turns and 
estimated as-deployed lengths at various points in the winding, 
including the estimated effects of Dyneema creep during bakeout, and the 
temperature and tension differences between winding and deployed 
conditions. (The deployed tension is assumed to be ~1.3 newtons.) The 
M L" in meters and "T" (in core turns) for each segment start at zero 
when that segment starts deploying. 

Quadratic length formula for PR0SED19 winding of F-2 tether: 

Noncond: 10132 meters; 14733 turns; L = 0.72822 * T - 2.75E-6 * Sqr(T) 

C-COR: 4865 meters; 11037 turns; L = 0.64497 * T - 18.5E-6 * Sqr(T) 

Insul: 215 meters; 1108 turns; L = 0.2339 * T - 36E-6 * Sqr(T). 

The F-2 non-conduct ive tether is 5 turns & 92 meters longer than F-l. 

The F-2 wire is 138 turns shorter but 12 meters longer than F-l. 

Most of these differences are due to the increased bulk of the longer 
overbraid on F-2, and the fact that we needed about the same wire length 
but the same # of Dyneema turns. 


Meters 

Turn# 

0 

0 

20 

28 

5562 

7872 

10132 

14733 

12705 

19328 

14334 

23384 

14997 

25770 

15012 

25835 

15212 

26878 


Transition 

Outer Kevlar tie-down 

End of Kevlar leader (start of Dyneema) 
Near the middle of the Dyneema 
Start of wire (end of Dyneema) 

Near the middle of the wire 
Parallel/criss-cross winding transition 
Start of overbraided insulated wire 
Brake turn (second of two to deploy) 
Termination pin at base of core 


F-2 cross-straps and key transitions (m) from inner termination pin: 


Meters 

CoreTurn 

Feature 

0 

0 

Termination pin at base of deployer 

200 

1043 

Brake enable sensor hole 

215 

1108 

End of overbraided insulated length. 

257 

1287 

1 

(inboard cross-strap; rest are spaced 321m apart) 

578 

2507 

2 


900 

3559 

3 


1221 

4 4 84 

4 


1542 

5331 

5 


1863 

6119 

6 


2184 

6853 

7 


2505 

7547 

8 


2826 

8206 

9 


3147 

8835 

10 


3469 

9438 

11 


3790 

10015 

12 


4111 

10572 

13 


4432 

11110 

14 


4753 

11631 

15 


5074 

12136 

16 

(outboard cross-strap) 

5080 

12145 

End 

of wire (7 wires tuck into core, 1 TT apart) 
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F-3 flight tether spooling data 

End-to-End Description of Flight Winding of F-3 (PROSED21) 

Joe Carroll 

Tether Applications, Inc. 

Winding dates: March 18-19, 2002 
Post-wind handling finished March 21, 2002 


Notes on winding: 


F-3 quadratic length formulas (winding 
Non-conduct : 10167 meters; 14730 turns; 
C-COR wire: 4870 meters; 11042 turns; 

Insul wire: 214 meters; 1103 turns; 


PR0SED21 ) : 

L = 0 . 73029*T - 2.72E-6*Sqr 
L = 0 . 64 642*T - 18 . 6E-6*Sqr 
L = 0 . 2337*T - 36E-6*Sqr(T) 


(T) 

(T) 


For comparison, here are the formulas for F-2 and F-l: 


F-2 quadratic length formulas (winding PROSED19) : 
Non-conduct: 10132 meters; 14733 turns; L = 0.72822*T 

C-COR wire: 4865 meters; 11037 turns; L = 0.64497*T 

Insul wire: 215 meters; 1108 turns; L = 0.2339*T 


2 . 75E-6*Sqr (T) 
18 . 5E-6*Sqr (T) 
36E-6*Sqr (T) . 


F-l quadratic 
Non-conduct : 
C-COR wire: 
Insul wire: 


length formulas (for first winding, PROSED15 ) : 

10040 meters; 14728 turns; L = 0.7215*T - 2 . 70E-6*Sqr (T) 

4894 meters; 11350 turns; L = 0.6378*T - 18 . 2E-6*Sqr (T) 

174 meters; 933 turns; L = 0.2201*T - 36E-6*Sqr (T) . 


Tether mass/length: 
Kevlar: 2. lOg/m 

Dyneema: 0.150g/m 

CCOR : 1 . 92g/m 

Inspire: 3.25g/m 


Transition locations for F-3 


As wound: 


MetW 

Core 

30917 

26875 

30878 

26849 

19596 

19008 

10278 

12145 

5072 

7556 

1777 

3497 

430 

1103 

402 

1044 

0 

0 


As deployed: 

Meters 

Turn# 

0 

0 

19 

26 

5577 

7867 

10167 

14730 

12741 

19319 

14371 

23378 

15037 

25772 

15051 

25831 

15251 

26875 


(PROSED21) : 

Transition 

Outer Kevlar tie-down 

End of Kevlar leader (start of Dyneema) 
Stop mid- Dyneema (0.2% longer in orbit) 
Start of wire (end of Dyneema) 

Stop mid-wire (wire=0.07% shorter) 
Parallel /criss-cross transition 
End of overbraid wrap-splice 
Brake turn (second of two to deploy) 
Termination pin at base of core 


Turn differences 
F-3: 


F-2 

F-l#l 

F-l# 3 

+ 3 

-2 

+ 76 

+ 6 

-233 

-154 

-2 

+ 306 

+ 453 

+ 4 

+ 102 

+ 228 

+ 3 

+ 136 

+ 272 


of F-2 (PR0SED1 9) and F-l ( PROSED15 & PROSED18) from 


Transition 

Start of wire (end of Dyneema) 
Parallel/criss-cross transition 
End of overbraid wrap-splice 
Brake turn (second of two to deploy) 
Termination pin at base of core 
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Basis of length estimates: . 

Estimated as-deployed length per turn of TAI metwheel (491 . 066mm/turn 
w/o tether) 496.5 mm for insulated overbraided tether & 20m Kevlar 
leader (same as as-wound) 494.5 mm for C-COR-coated conductive tether 
(this is 0.07% shorter than as-wound) 492.6 mm for flat 13x100 Dyneema 
braid (this is 0.2% longer than as-wound) 

Winding temperature range: -21-23C 

Winding tension (lengths were measured at these tensions) : 

13-16 newtons for parallel winding 

12, 10.5, 9, 7.5, and 10 newtons for CC patterns #1-5. 

Explanation of adjustments from as-wound to as-deployed conditions: 
Assume Dyneema creep during bakeout reduces its tension to as-deployed 

tension. , „ . 

Assume Kevlar has equilibrium post-deployment tension of -1.3 newtons 
Modulus of long segments (C-COR and Dyneema) is -12000-15000 newtons. 
Assume Dyneema and CCOR are 225K and 300K in orbit; CTEs are -28 PP m/K 
and 0 P pm/K . Thus Dyneema should be -0.2% longer; CCOR -0.07% shorter, 
deployed vs. wound. 



The differences in turns of the F-2 from the F-3 are less than 6 turns which makes them 
almost perfectly alike to the control law. The reference profile (Ref. #78) works equally 
well for the F-2 and F-3 flight tethers. This reference profile, shown in Figure 17, has 
several good features as follows: (a) the exit velocity changes smoothly during the Dyneema 
deployment and it does not deep down to low values; (b) the brake increases monotonously 
from zero to its maximum value; and (c) the libration angle is very close to zero at the end of 
deployment (for reference conditions). This reference profile when combined with the 
control law provides a final libration amplitude which is quite insensitive to variation of the 
Dyneema minimum tension and other friction parameters as it will be shown later on. The 
reference table for Ref. #78 is tabulated in Appendix A. 
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Reference Profile #78 



Figure 17 Control law Reference Profile #78 


3.4 Flight control law parameters 

The control law not only provides the logic for controlling deployment under perfect 
functionality of the sensors. It also has a set of backup modes to circumvent or mitigate the 
effects of possible failures or malfunctions of the deployer sensors. The sensors in 
question are the two turn counters (Counter- A and Counter-B) which measure the number 
of tether turns that have been deployed. In standard operation, the output of the two turn 
counters is combined to produce a software turn count Counter-C which is unaffected by 
double counts. In case of failure of one of the turn counter, Counter-C is bypassed and the 
still functioning turn counter is utilized by the sensor logic. However, when only one turn 
counter is working there is an over count of the number of turns that has to be accounted for 
by the control law when transitioning from one tether segment to the next. Also, if both turn 
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counters fail, the control law logic sets the feedback to zero and utilizes exclusively the 
reference profile to drive deployment (this is one additional advantage of having the feed- 
forward profile in the control law). Similarly, there are two levels of backup for the 
actuation of the slow down maneuver that starts at the beginning of the insulated tether 
portion (last 215 m of tether). The slow down maneuver is activated, in standard operation, 
by the turn counter. In case of the two turn counters failing, there are two level of backup as 
follows: the first one is provided by the brake enabler switch which senses the beginning of 
the insulated wire and the second is provided by time. 

The control parameters define the standard operation of the control law and also the 
activation of the backup modes. The parameters controlling the standard operation are taken 
from the SEDS-II control set with a few adjustments aimed at improving the robustness of 
the control law to variations of the friction parameters. The control parameters that define 
the control modes corresponding to transitions from one tether segment to the next and the 
parameters that define the backup modes have been added to the original control set of the 
SEDS-II control law. The complete set of the updated control parameters for Reference 
#78 are shown in Table 1 . 
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Table 1 Updated control parameters for Reference #78 


Values of control parameters for ProSEDS Ref#78 (F2 & F3 
tethers). Issued on April 5, 2002; revised on July 23 , 2002. 


CONTROL PARAMETERS 


No. 

PARAMETER 

VALUE 

(Units) 

1. 

C 

0.125 


2. 

GTC 

0.002 

(1/turn) 

3. 

DZTC 

5 

( turn ) 

4. 

TCELIM 

3000 

(turn) 

5. 

GTCR 

0.4 

( s/turn) 

6. 

DZTCR 

0.1 

(turn/s) 

7. 

TCRELIM 

5 

(turn/s) 

8. 

WAILP 

3 


9. 

TBD s 

65535 

( sec ) 

10. 

BIAS 

0 

(turn) 

11. 

WACLP 

6.5 

( turn ) 

12. 

TCBS 

18000 

(turn) 

13. 

A1 

0.724 


14. 

A2 

2.82E- 

■6 

15. 

STOPDEPLOY 

65535 

(sec) 

16. 

TCDUTY 

13560 

(turn) 

17. 

TURNBRAKE0 

14265 

(turn) 

CO 

• 

WIREBRAKE 

0.5 

(turn) 

19. 

RAMPUP 

25835 

(turn) 


(REF. #78) 

Type 

Filter coefficient 

TurnCount Gain 

TurnCount Deadzone 

Max. TurnCount Error 

TurnCountRate Gain 

TurnCountRate Deadzone 

Max. TurnCountRate Error 

Wraplncrement UpperLimit 

Time after which BIAS is 
applied 

BrakePost Bias 

WrapAdjustment UpperLimit 

Turns Count Brake Stop 
(SEDS-II law) 

Coeff_l in Variable Gains 

Coeff 2 in Variable Gains 

Time for brake ramping up 
to max. value at end of 
deployment (SEDS-II law) 

End of 50% duty cycle 

RampDown brake to 
WIREBRAKE 

BrakeTurns during CCOR 
deployment 

Start of slowdown 
procedure 
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20. QUITLAWBACKUP 14400 (turn) 


21. BRSD 

1.25 

(turn) 

22. TBD( 15 ) 

20 

(sec) 

23. TIMECFAIL 

120 

(sec) 

24. TIMEDUTY 

6080* (sec) 

25. TIMEQUERY 

6832* (sec) 


26. TIMERAMPNOBES 6894* (sec) 


Ramp down brake in case 
of CounterA or B failure 

Max brake turns during 
slow down 

Time to rampup brake from 
WIREBRAKE to BRSD 

Time of no update of 
CounterC to declare the 
CounterC failed 

Time-based equivalent of 
TCDUTY 

The software starts 
interrogating the BES if 
CounterC has failed 

Time-based start of 
slowdown procedure if the 
BES was declared failed 


♦These values of parameters assume a SMET = 2640 sec at 
endmass separation (parameters revised on July 23, 2002) 


Notes : 

BES = brake enable switch 

SW = flight software 

SMET = secondary mission elapsed time 


3.5 Sensitivity of control law to friction parameters 

Deployment trajectory is mostly influenced by the tether tension profile during the early 
part of deployment (see Ref. 8 ). The physical explanation for this behavior is that small 
changes in the tether changes alter the dynamics the most when the external forces (i.e., the 
gravity gradient) are small, that is, at short tether lengths. 

The deployment test data show that the range of variation of most of the faction 
parameters are smaller than ±20%. The minimum tether tensions of the Dyneema and the 
CCOR wire are the exception. The value of the Dyneema minimum tension in particular 
changes by a few folds across the many tether samples tested. The CCOR minimum 
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tension, besides being less variable than the Dyneema, has an almost negligible effect on the 
libration amplitude at the end of deployment. 

The minimum tension of the ProSEDS Dyneema tether (which dominates the final state 
at the end of deployment) has already been measured in deployment ground tests under 
different temperatures to vary between 5mN and 30 mN for all the tether samples tested 
with a smaller range of variation for the F-l tether. Consequently, the control law must 
provide a residual libration at the end of deployment of less than 20° (as specified by the 
mission requirements) within the measured range of variability of the minimum tension. 

The control law with Ref. #78 can tolerate without a significant decay in performance a 
value of the non-conductive tether minimum tension between 5 mN and 30 mN. For 30 mN 
< T 0 < 60 mN, the libration at end of deployment increases above 10° but still below 20°. 
For T 0 & 80 mN, the deployment stops at a distance of about 500 m because of excessive 
friction and without any role being played by the control law. The critical value of 80 mN 
for the minimum tension is determined by the ejection velocity which with the present 
ejection system is equal to 2.74 m/s. 

Figure 18-20 show details of the deployment dynamics for values of the minimum 
Dyneema tension for T 0 = 5 mN, 10 mN, and 20 mN (reference value), respectively. 
Figure 21 shows the effect on the control law and deployment dynamics of a noisy tether 
tension on the deployment dynamics (for T 0 = 20 mN). A ±50% white noise has been 
added to the tether tension. The figure shows that the control law is capable of handling the 
high noise level with ease. The brake actuation is still smooth and the deployment dynamics 
is negligibly affected by the noise. The final libration amplitude is unaffected by the white 
noise. 
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Ref#78, Tmin = 5mN/1 50mN, AV = 2.74 m/s, Brake 1 ,25/20s 



Figure 18 Deployment dynamics for Ref#78 and T 0 = 5 mN 
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Ref#78, Tmin = 10mN/150mN, AV = 2.74 m/s, Brake 1 ,25/20s 



Figure 19 Deployment dynamics for Ref#78 and T 0 = 10 mN 
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Ref#78, Tmin = 20mN/1 50mN, AV = 2.74 m/s, Brake 1 .25/20s 



Figure 20 Deployment dynamics for Ref#78 and T 0 = 20 mN (reference) 
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Ref#78, Tmin = 20mN/150mN, AV = 2.74 m/s, ±50% White Noise 



Figure 21 Deployment dynamics for T 0 = 20 mN with ±50% tension white noise 
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Local horizontal (km) 

Figure 22 Deployment trajectory of the endmass with respect to the Delta stage 

Figure 22 depicts the deployment trajectory of the endmass with respect to the Delta 
stage for the reference case of T 0 = 20 mN. The tick marks along the trajectory are every 6 
minutes and the CCOR wire starts deploying at about the point where the tick mark caption 
is pointing. Figure 23 shows the amplitude of the residual libration at the end of 
deployment vs. the minimum tension T 0 of the Dyneema tether for Ref. #78. The Figure 
also shows the libration amplitude without deployment control and the representative points 
of the SEDS-I and SEDS-II missions. The final libration amplitude is sensitive to the 
Dyneema tether T 0 and it is quite insensitive to the value of the wire minimum tension T wjre . 
Values of T wire of 50-300 mN have been explored with very good deployment dynamics. 
Values as high as 500 mN are tolerable for the minimum tension of the wire. 
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Figure 23 Final libration amplitude vs. T 0 with control (Ref. #78) and without control 


3.6 Concluding remarks 

The deployment control law in combination with Reference Profile #78 is robust with 
respect to changes of the most influential and uncertain friction parameters. A number of 
backup modes have been built into the control logic to circumvent (or mitigate) the effects of 
the possible failure of one (or two) tum counters. The control law meets the requirement of 
a final libration smaller than 20° over the measured values of the friction parameters with a 
positive safety margin. 
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4.0 Analysis/Estimation of deployment flight data 

4.1 Introduction 

The performance of the deployment control law in terms of final libration (and other 
dynamic variables) can be assessed from the single data stream of the turn counter. In fact, 
from the (in-plane) libration equation it is possible to see that in the absence of strong 
external forces (other than gravity gradient forces) the libration state of the system can be 
extracted from the knowledge of the tether length and speed profiles. Since deployment 
takes place almost exclusively on the orbital plane, knowledge of the in-plane angle define 
the libration state of the system. 

It is possible to construct an observer of the system libration based on the deployed 
tether length and exit speed. In the actual case, the only available information is the turn 
counter which provides the number of deployed tether turns. The tether turns must be first 
converted into tether length by utilizing the spooling relationship. The length data must be 
filtered with a low-pass filter before computing the exit speed because any small jump in 
tether length over a short time period will result in a big jump of tether rate when the speed 
is computed by numerical derivation. 

Finally, the smoothed length and speed can be used as input to the deployment observer 
in order to compute tether libration and consequently reconstruct the deployment trajectory 
of the endmass and estimate the libration amplitude after the end of deployment when the 
attitude of the endmass has stabilized. 

An alternative way to the deployment state estimation is to utilize position estimates 
from the GPS receivers on board the endmass and the Delta and derive the relative 
dynamics through differences of simultaneous position data. The GPS-method, however, is 
predicated on having the two signals available from the very beginning of deployment when 
the endmass is tumbling and prone to loosing the signal lock. The GPS method, which will 
be implemented later on, could be utilized more reliably to backup the estimate of the 
libration amplitude at the end of deployment. 
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4.2 Numerical results 

In the following we will show results from the deployment observer (i.e., based on the 
turn counter estimation method). We take a noisy string of turn counter data from the 
deployment simulation shown in Fig. 21 (with T 0 = 20 mN and heavy white noise) and we 
prove that the deployment trajectory can be reconstructed and the amplitude of the final 
libration accurately evaluated by means of the process described previously. 



LH (km) 

Figure 24 Deployment trajectory of endmass: estimated from noisy turn counter data 
(solid line); and original data (dotted line) 
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Figure 24 shows the deployment trajectory of the endmass with respect to the Delta 
stage estimated from noisy tum counter data (for T 0 = 20 mN) and the original trajectory 
computed from the simulated data shown in Fig. 21. Clearly, the estimation software is 
capable of extracting the deployment trajectory from noisy tum counter data and of 
estimating accurately the libration amplitude at the end of deployment. This estimation 
software was also installed on the computer at NASA MSFC and it was run successfully 
twice during the mission operation simulations at the EDAC. 


4.3 Concluding remarks 

The estimation technique based on the data string provided by the tum counter is able to 
reconstruct the deployment trajectory and provide an accurate estimate of the amplitude of 
the libration at the end of deployment starting from noisy data. The advantage of this 
method is that it is based on one string of data that is readily available throughout 
deployment. The GPS data from both the Delta stage and the endmass can also be utilized 
to check the estimate of the libration after the end of deployment. 
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Appendix A 

Reference Table for ProSEDS deployment control law (Ref#78, flight). For completeness, 
this file includes data of tether length and tether exit velocity 


Date: April 4, 2002 
Author: E.C. Lorenzini (SAO) 

Support File for ProSEDS Reference #78. 


Time 

Length 

Speed 

Is) 

(m) 

(m/s) 

0 

0.5 

2.74 

8 

22.378 

2.73 

16 

44 . 173 

2.719 

24 

65.889 

2.71 

32 

87.53 

2.701 

40 

109. 101 

2.692 

48 

130.605 

2.684 

56 

152.044 

2.676 

64 

173.424 

2.669 

72 

194.749 

2.662 

80 

216.021 

2.656 

88 

237,245 

2.65 

96 

258.424 

2.645 

104 

279.561 

2.64 

112 

300.662 

2.635 

120 

321.729 

2 . 631 

128 

342.765 

2.628 

136 

363.775 

2.625 

144 

384.761 

2.622 

152 

405.728 

2.62 

160 

426.679 

2.618 

168 

447.617 

2 . 617 

176 

468.545 

2.616 

184 

489.467 

2.615 

192 

510.386 

2.615 

200 

531.306 

2.615 

208 

552.229 

2.616 

216 

573.159 

2.617 

224 

594 .098 

2.618 

232 

615.05 

2.62 

240 

636.019 

2.622 

248 

657.005 

2 . 625 

256 

678 . 014 

2.628 

264 

699.046 

2.631 

272 

720.106 

2.634 

280 

741 . 196 

2 . 638 

288 

762.318 

2.642 

296 

783.475 

2.647 

304 

804 . 67 

2.652 

312 

825.904 

2.657 


Turn 

Turnrate 

Brake 

(turn) 

(turn/s) 

(turn) 

1 

3.752 

0 

31 

3.739 

0 

61 

3.726 

0 

90 

3.714 

0 

120 

3.703 

0 

150 

3.692 

0 

179 

3.681 

0 

208 

3.672 

0 

238 

3.663 

0 

267 

3.654 

0 

296 

3.647 

0 

325 

3. 64 

0 

354 

3.633 

0 

384 

3.627 

0 

413 

3.622 

0 

441 

3.618 

0 

470 

3.614 

0 

499 

3.61 

0 

528 

3.608 

0 

557 

3.605 

0 

586 

3.604 

0 

615 

3.603 

0 

643 

3.602 

0 

672 

3.602 

0 

701 

3.603 

0 

730 

3.604 

0 

759 

3.606 

0 

788 

3.609 

0 

817 

3.611 

0 

845 

3.615 

0 

874 

3.619 

0 

903 

3.623 

0 

932 

3.628 

0 

961 

3.633 

0 

990 

3.639 

0 

1020 

3.645 

0 

1049 

3.652 

0 

1078 

3.659 

0 

1107 

3.667 

0 

1137 

3.675 

0 
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320 

847.182 

2.662 

1166 

3.684 

0 

328 

868.504 

2.668 

1196 

3.692 

0 

336 

889.873 

2.674 

1225 

3.702 

0 

344 

911.291 

2.681 

1255 

3.711 

0 

352 

932.76 

2.687 

1285 

3.721 

0 

360 

954.284 

2.694 

1314 

3.732 

0 

368 

975.863 

2.701 

1344 

3.743 

0 

376 

997 .499 

2.708 

1374 

3.754 

0 

384 

1019.194 

2.716 

1404 

3.765 

0 

392 

1040.951 

2.724 

1435 

3.777 

0 

400 

1062.771 

2.731 

1465 

3.789 

0 

408 

1084 . 655 

2.74 

1495 

3.801 

0 

416 

1106.606 

2.748 

1526 

3.814 

0 

424 

1128.625 

2.757 

1556 

3.827 

0 

432 

1150.713 

2.765 

1587 

3.84 

0 

440 

1172.872 

2.774 

1618 

3.853 

0 

448 

1195.103 

2.783 

1649 

3,867 

0 

456 

1217.408 

2.793 

1680 

3.881 

0 

464 

1239.787 

2.802 

1711 

3.895 

0 

472 

1262.242 

2.812 

1742 

3.909 

0 

480 

1284.775 

2.821 

1773 

3.924 

0 

488 

1307 . 386 

2.831 

1805 

3.939 

0 

496 

1330.076 

2.841 

1836 

3.954 

0 

504 

1352 .846 

2.851 

1868 

3.969 

0 

512 

1375.697 

2.862 

1900 

3.984 

0 

520 

1398.631 

2.872 

1932 

4 

0 

528 

1421 . 646 

2.882 

1964 

4.015 

0 

536 

1444 .746 

2.893 

1996 

4.031 

0 

544 

1467 . 929 

2.903 

2028 

4 . 047 

0 

552 

1491 . 196 

2.914 

2061 

4.062 

0 

560 

1514 . 549 

2.924 

2093 

4.079 

0 

568 

1537 . 988 

2 . 935 

2126 

4.095 

0 

576 

1561 . 513 

2.946 

2159 

4 .111 

0 

584 

1585.124 

2 . 957 

2192 

4 . 127 

0 

592 

1608 . 822 

2.968 

2225 

4 . 143 

0 

600 

1632.607 

2.979 

2258 

4 . 16 

0 

608 

1656.48 

2.99 

2291 

4 . 176 

0 

616 

1680.44 

3 

2325 

4 .193 

0 

624 

1704 . 488 

3.011 

2358 

4 . 209 

0 

632 

1728 . 623 

3.022 

2392 

4.226 

0 

640 

1752 . 846 

3.033 

2426 

4 .242 

0 

648 

1777.157 

3.044 

2460 

4 .259 

0 

656 

1801.556 

3.055 

2494 

4.275 

0 

664 

1826.042 

3.066 

2528 

4.292 

0 

672 

1850. 615 

3.077 

2563 

4.309 

0 

680 

1875.276 

3.088 

2597 

4.325 

0 

688 

1900.024 

3.099 

2632 

4.342 

0 

696 

1924 . 859 

3.11 

2667 

4 . 358 

0 

704 

1949.78 

3. 121 

2702 

4.374 

0 

712 

1974 .787 

3.131 

2737 

4.391 

0 

720 

1999.88 

3.142 

2772 

4 .407 

0 

728 

2025.058 

3.153 

2807 

4 .423 

0 

736 

2050.321 

3.163 

2843 

4 .44 

0 

744 

2075.669 

3.174 

2878 

4 .456 

0 

752 

2101. 1 

3.184 

2914 

4.472 

0 

760 

2126.615 

3.194 

2950 

4 .487 

0 

768 

2152.212 

3.205 

2986 

4.503 

0 

776 

2177 .891 

3.215 

3022 

4.519 

0 
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784 

2203.652 

3.225 

792 

2229.493 

3,235 

800 

2255.415 

3.245 

808 

2281.416 

3.255 

816 

2307.495 

3.265 

824 

2333.652 

3.274 

832 

2359.886 

3.284 

840 

2386.197 

3.293 

848 

2412.582 

3.303 

856 

2439.042 

3.312 

864 

2465.575 

3.321 

872 

2492.181 

3.33 

880 

2518.859 

3.339 

888 

2545.607 

3.348 

896 

2572.425 

3.357 

904 

2599.312 

3.365 

912 

2626.267 

3.374 

920 

2653.287 

3.382 

928 

2680.372 

3.39 

936 

2707 . 52 

3.397 

944 

2734.729 

3.405 

952 

2761.997 

3.412 

960 

2789.323 

3.419 

968 

2816.704 

3.426 

976 

2844 . 14 

3.433 

984 

2871 . 628 

3.439 

992 

2899.166 

3.445 

1000 

2926.752 

3.451 

1008 

2954.384 

3.457 

1016 

2982.061 

3.4 62 

1024 

3009.78 

3.4 67 

1032 

3037 . 54 

3.472 

1040 

3065.338 

3.477 

1048 

3093.173 

3.481 

1056 

3121.042 

3.486 

1064 

3148 . 943 

3.49 

1072 

3176.875 

3.493 

1080 

3204 .836 

3.497 

1088 

3232.822 

3.5 

1096 

3260.833 

3.503 

1104 

3288.866 

3.505 

1112 

3316.919 

3.508 

1120 

3344 . 99 

3.51 

1128 

3373.076 

3.512 

1136 

3401.177 

3.513 

1144 

3429.29 

3.515 

1152 

3457.412 

3.516 

1160 

3485.542 

3.517 

1168 

3513.677 

3.517 

1176 

3541.816 

3.517 

1184 

3569.955 

3.517 

1192 

3598.094 

3.517 

1200 

3626.23 

3.517 

1208 

3654.361 

3.516 

1216 

3682.485 

3.515 

1224 

3710.6 

3.514 

1232 

3738.704 

3.512 

1240 

3766.794 

3.51 


3058 

4.535 

0 

3094 

4 . 55 

0 

3131 

4.566 

0 

3168 

4.581 

0 

3204 

4.596 

0 

3241 

4 . 611 

0 

3278 

4 . 626 

0 

3315 

4 . 641 

0 

3352 

4 , 655 

0 

3390 

4 . 67 

0 

3427 

4 . 684 

0 

3465 

4.698 

0 

3502 

4.712 

0 

3540 

4.726 

0 

3578 

4 . 74 

0 

3616 

4.754 

0.004 

3654 

4.767 

0.028 

3692 

4.78 

0.048 

3730 

4.793 

0.067 

3769 

4.805 

0.086 

3807 

4.817 

0.104 

384 6 

4 . 829 

0.123 

3884 

4.841 

0.142 

3923 

4.852 

0.16 

3962 

4.863 

0.179 

4001 

4.874 

0.197 

4040 

4.884 

0.215 

4079 

4 .894 

0.233 

4118 

4 . 904 

0.251 

4158 

4.913 

0.268 

4197 

4 . 922 

0.285 

4236 

4 . 93 

0.303 

4276 

4 . 939 

0.32 

4315 

4.946 

0.337 

4355 

4 . 954 

0.353 

4395 

4 . 961 

0.37 

4434 

4 . 968 

0.387 

4474 

4 . 975 

0.403 

4514 

4 . 981 

0.42 

4554 

4 . 986 

0 .436 

4594 

4 . 992 

0.452 

4634 

4 . 997 

0.469 

4674 

5.001 

0.485 

4714 

5.006 

0.501 

4754 

5.01 

0.517 

4794 

5.013 

0.533 

4834 

5.016 

0.549 

4874 

5.019 

0.565 

4914 

5.022 

0.581 

4955 

5.024 

0.597 

4995 

5.025 

0.612 

5035 

5.027 

0.628 

5075 

5.028 

0 . 644 

5115 

5.028 

0.659 

5156 

5.028 

0.675 

5196 

5.028 

0.69 

5236 

5.028 

0.706 

5276 

5.027 

0.721 
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1248 

3794.868 

3.508 

1256 

3822.926 

3.506 

1264 

3850.964 

3.503 

1272 

3878 . 981 

3.501 

1280 

3906.974 

3.498 

1288 

3934 . 941 

3.494 

1296 

3962.881 

3.491 

1304 

3990.791 

3.487 

1312 

4018.67 

3.483 

1320 

4046.516 

3.478 

1328 

4074 . 326 

3.474 

1336 

4102.099 

3.469 

1344 

4129.833 

3.464 

1352 

4157,525 

3.459 

1360 

4185.175 

3.453 

1368 

4212.779 

3.448 

1376 

4240.337 

3.442 

1384 

4267.847 

3.436 

1392 

4295.307 

3.429 

1400 

4322,714 

3.423 

1408 

4350.068 

3.416 

1416 

4377 . 366 

3.409 

1424 

4404 . 607 

3.401 

1432 

4431.789 

3.394 

1440 

4458 . 911 

3.386 

1448 

4485.971 

3.378 

1456 

4512.967 

3.371 

1464 

4539.898 

3.362 

1472 

4566.763 

3.354 

1480 

4593.559 

3.345 

1488 

4620.285 

3.336 

1496 

4646.941 

3.327 

1504 

4673.524 

3.318 

1512 

4700.033 

3.309 

1520 

4726.468 

3.299 

1528 

4752.825 

3.29 

1536 

4779.106 

3.28 

1544 

4805.307 

3.27 

1552 

4831 .429 

3.26 

1560 

4857.469 

3.25 

1568 

4883.428 

3.24 

1576 

4909.303 

3.229 

1584 

4935.094 

3.219 

1592 

4960.799 

3.208 

1600 

4986.419 

3.197 

1608 

5011.951 

3.186 

1616 

5037.396 

3. 175 

1624 

5062.753 

3.164 

1632 

5088.02 

3.153 

1640 

5113.197 

3.141 

1648 

5138.283 

3.13 

1656 

5163.279 

3.119 

1664 

5188 . 182 

3.107 

1672 

5212.994 

3.096 

1680 

5237.713 

3.084 

1688 

5262.339 

3.072 

1696 

5286.871 

3.061 

1704 

5311.31 

3.049 


5316 

5.025 

0.737 

5357 

5.024 

0.752 

5397 

5.022 

0.767 

5437 

5.019 

0.782 

5477 

5.017 

0.798 

5517 

5.013 

0.813 

5557 

5.01 

0.828 

5597 

5.006 

0.843 

5637 

5.002 

0.858 

5677 

4 . 997 

0.873 

5717 

4 , 992 

0.888 

5757 

4 . 987 

0.903 

5797 

4 . 981 

0.917 

5837 

4 . 975 

0.932 

5877 

4.969 

0.947 

5917 

4.962 

0.962 

5956 

4 . 956 

0.976 

5996 

4 . 948 

0.991 

6035 

4 . 94 

1.005 

6075 

4 . 933 

1.02 

6114 

4 . 924 

1.034 

6154 

4 . 916 

1.049 

6193 

4 . 907 

1.063 

6232 

4.897 

1.077 

6271 

4 . 888 

1.091 

6310 

4 . 878 

1.105 

6349 

4 .868 

1,119 

6388 

4.858 

1 . 134 

6427 

4.847 

1.148 

64 66 

4.836 

1.161 

6504 

4.825 

1 . 175 

6543 

4.813 

1.189 

6581 

4.802 

1.203 

6620 

4.79 

1.217 

6658 

4.777 

1.23 

6696 

4.765 

1.244 

6734 

4.752 

1.257 

6772 

4.739 

1.271 

6810 

4.726 

1.284 

6848 

4.712 

1.298 

6886 

4 . 699 

1.311 

6923 

4 . 685 

1.324 

6961 

4 . 671 

1.337 

6998 

4 . 657 

1.35 

7035 

4.643 

1.363 

7072 

4 . 628 

1.376 

7109 

4 . 613 

1.389 

7146 

4.599 

1.402 

7183 

4.584 

1.415 

7219 

4 . 569 

1 . 427 

7256 

4 . 553 

1.44 

7292 

4.538 

1 .452 

7328 

4.523 

1.465 

7364 

4.507 

1.477 

7400 

4.491 

1.489 

7436 

4 .476 

1.502 

7472 

4 .46 

1.514 

7508 

4 .444 

1.526 
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1712 

5335.655 

3.037 

1720 

5359.906 

3.026 

1728 

5384 . 063 

3.014 

1736 

5408.126 

3.002 

1744 

5432 . 095 

2.99 

1752 

5455.969 

2.978 

1760 

5479.75 

2.967 

1768 

5503.436 

2.955 

1776 

5527.029 

2.943 

1784 

5550.528 

2.932 

1792 

5573.935 

2.92 

1800 

5597.248 

2,908 

1808 

5620.461 

2.895 

1816 

5643.564 

2.881 

1824 

5666.556 

2.867 

1832 

5689.439 

2.854 

1840 

5712.214 

2.84 

1848 

5734.88 

2.826 

1856 

5757.439 

2.813 

1864 

5779.891 

2.8 

1872 

5802.237 

2.787 

1880 

5824 .479 

2.774 

1888 

5846.616 

2.761 

1896 

5868.651 

2.748 

1904 

5890.583 

2.735 

1912 

5912.415 

2.723 

1920 

5934 . 148 

2.71 

1928 

5955.782 

2.698 

1936 

5977 . 32 

2.686 

1944 

5998.762 

2.674 

1952 

6020.11 

2.663 

1960 

6041.365 

2. 651 

1968 

6062 . 53 

2.64 

1976 

6083.605 

2.629 

1984 

6104 . 593 

2.618 

1992 

6125.495 

2.608 

2000 

6146.313 

2.597 

2008 

6167.049 

2.587 

2016 

6187.704 

2.577 

2024 

6208.282 

2.567 

2032 

6228.783 

2.558 

2040 

6249.21 

2.549 

2048 

6269.566 

2.54 

2056 

6289.852 

2.532 

2064 

6310.07 

2.523 

2072 

6330.223 

2.515 

2080 

6350.315 

2.508 

2088 

6370.345 

2.5 

2096 

6390.319 

2.493 

2104 

6410.236 

2.486 

2112 

6430.102 

2.48 

2120 

6449. 918 

2.474 

2128 

6469.687 

2.468 

2136 

6489.411 

2.463 

2144 

6509.094 

2.458 

2152 

6528.738 

2.453 

2160 

6548.347 

2.449 

2168 

6567.923 

2.445 


7543 

4.428 

1.538 

7579 

4.412 

1.55 

7614 

4.396 

1.562 

7649 

4.38 

1.574 

7684 

4.364 

1.586 

7719 

4 . 348 

1.597 

7753 

4.332 

1.609 

7788 

4.316 

1.621 

7822 

4.301 

1.632 

7857 

4.285 

1 . 644 

7891 

4.269 

1.656 

7925 

4.253 

1.667 

7959 

4.234 

1.678 

7993 

4.215 

1.69 

8027 

4 .196 

1.701 

8060 

4 . 177 

1.713 

8093 

4 . 158 

1.724 

8127 

4 . 14 

1.735 

8160 

4 . 121 

1.746 

8192 

4.103 

1.758 

8225 

4.084 

1.769 

8258 

4.066 

1.78 

8290 

4 . 048 

1.791 

8323 

4.031 

1.802 

8355 

4.013 

1.813 

8387 

3.996 

1.824 

8419 

3.979 

1.835 

8450 

3.962 

1.846 

8482 

3.945 

1.857 

8514 

3.928 

1 .868 

8545 

3.912 

1 . 879 

8576 

3.896 

1.89 

8607 

3.881 

1.901 

8638 

3.865 

1.912 

8669 

3.85 

1.923 

8700 

3.836 

1 . 934 

8731 

3.821 

1 . 945 

8761 

3.807 

1.956 

8791 

3.794 

1.967 

8822 

3.78 

1.978 

8852 

3.767 

1 . 988 

8882 

3.755 

1 . 999 

8912 

3.743 

2.01 

8942 

3.731 

2.021 

8972 

3.719 

2.032 

9001 

3.709 

2.043 

9031 

3.698 

2.054 

9061 

3.688 

2.065 

9090 

3.678 

2.076 

9119 

3.669 

2.087 

9149 

3.661 

2.098 

9178 

3.653 

2.109 

9207 

3.645 

2.12 

9236 

3.638 

2 . 131 

9265 

3.631 

2.142 

9294 

3.625 

2.153 

9323 

3.62 

2.164 

9352 

3.615 

2.175 
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2176 

6587.469 

2.441 

2184 

6606.988 

2.438 

2192 

6626.484 

2.436 

2200 

6645 . 96 

2.433 

2208 

6665.419 

2.431 

2216 

6684 .863 

2 .43 

2224 

6704 .296 

2.429 

2232 

6723.723 

2.428 

2240 

6743. 145 

2.428 

2248 

6762.567 

2.428 

2256 

6781 . 991 

2.428 

2264 

6801 .421 

2.429 

2272 

6820.861 

2.431 

2280 

6840.315 

2.433 

2288 

6859.785 

2 .435 

2296 

6879.275 

2.438 

2304 

6898.789 

2.441 

2312 

6918 . 331 

2.445 

2320 

6937 . 903 

2.449 

2328 

6957 . 51 

2 .453 

2336 

6977 . 153 

2.458 

2344 

6996.836 

2.463 

2352 

7016.562 

2.468 

2360 

7036.332 

2 .474 

2368 

7056.149 

2.48 

2376 

7076.014 

2 .486 

2384 

7095 . 93 

2.493 

2392 

7115 .898 

2.499 

2400 

7135.919 

2 . 506 

2408 

7155 . 993 

2 .513 

2416 

7176. 123 

2.52 

2424 

7196,308 

2.527 

2432 

7216.549 

2.534 

2440 

7236.845 

2.54 

2448 

7257 .197 

2.547 

2456 

7277 . 604 

2.554 

2464 

7298 . 066 

2.561 

2472 

7318.581 

2.568 

2480 

7339.148 

2.574 

2488 

7359.765 

2.58 

2496 

7380.432 

2.586 

2504 

7401 . 145 

2.592 

2512 

7421 . 902 

2.597 

2520 

7442.701 

2.602 

2528 

7463.539 

2.607 

2536 

7484 .412 

2.611 

2544 

7505.318 

2.615 

2552 

7526.251 

2.618 

2560 

7547.209 

2.621 

2568 

7568.186 

2 . 623 

2576 

7589.179 

2.625 

2584 

7610.183 

2.626 

2592 

7631 . 19 

2.626 

2600 

7652 . 198 

2.626 

2608 

7673.201 

2.625 

2616 

7694 . 192 

2.623 

2624 

7715.165 

2.62 

2632 

7736.115 

2.617 


9381 

3.61 

2.186 

9410 

3.606 

2.197 

9439 

3.603 

2.209 

9468 

3.601 

2.22 

9497 

3.598 

2.231 

9525 

3.597 

2.242 

9554 

3.596 

2.254 

9583 

3.596 

2.265 

9612 

3.596 

2.276 

9640 

3.597 

2.288 

9669 

3.599 

2.299 

9698 

3.601 

2.311 

9727 

3.604 

2.322 

9756 

3.607 

2.334 

9785 

3.612 

2.345 

9813 

3.616 

2.357 

9842 

3.622 

2.369 

9871 

3.628 

2.381 

9900 

3.635 

2.394 

9930 

3.642 

2.408 

9959 

3.65 

2 .423 

9988 

3.659 

2.439 

10017 

3.668 

2.455 

10047 

3.677 

2 .472 
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3336 9134.229 
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3448 9393.606 

3456 9414.015 
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